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ABSTRACT 

We present the mid-infrared colors of X-ray-detected AGN and explore mid-infrared selection cri- 
teria. Using a statistical matching technique, the likelihood ratio, over 900 IRAC counterparts were 
identified with a new MUSYC X-ray source catalog that includes ^ 1000 published X-ray sources in 
the Chandra Deep Field-South and Extended Chandra Deep Field-South. Most X-ray-selected AGN 
have IRAC spectral shapes consistent with power-law slopes, fy cx t"", and display a wide range of 
colors, —2 < a < 2. Although X-ray sources typically fit to redder (more negative a) power-laws 
than non- X-ray detected galaxies, more than 50% do have flat or blue (galaxy-like) spectral shapes 
in the observed 3-8 /im band. Only a quarter of the X-ray selected AGN detected at 24 /im are well 
fit by featureless red power laws in the observed 3.6-24 /xm, likely the subset of our sample whose 
infrared spectra are dominated by emission from the central AGN region. Most IRAC color-selection 
criteria fail to identify the majority of X-ray-selected AGN, finding only the more luminous AGN, 
the majority of which have broad emission lines. In deep surveys, these color-selection criteria select 
10-20% of the entire galaxy population and miss many moderate luminosity AGN. 
Subject headings: (galaxies:) quasars: general — infrared: galaxies — X-rays: galaxies 



1. INTRODUCTION 

Theoretical arguments and observational evidence in- 
dicate that Active Galactic Nuclei ( AGN) may play an 
important role in galaxy evolution ( iFabian et al. | [20031. 
2006; Croton et al. 200 a iKriek et al . 2006; Siiacki^et^ 
20071: ISchawinski et al.H2007D The studv of this mutual 
feedback requires a complete sample of AGN, their host 
galaxies, and comparable inactive galaxies at the same 
epochs. To discover these populations, surveys with deep 
multi-wavelength coverage are necessary. 

Today, AGN are identified at many wavelengths, 
each of which selects a slightly different sample. 
For exampl e, radio surveys prefer entially find radio- 
loud AGN (jBrinkmann et al.ll2000| ). while optical and 
ultrav iolet-excess surveys find luminous unobscure d 
AGN (jRichards et al.ll2004[Mazzarella fc Balzano" 1986h . 
Because they detect all but the most heavily ob- 
scured AGN, hard X-ray surveys efficiently select 
relatively unbiased and complete samples of AGN 
(jBrandt fc Hasingeil 120051 ). even if X-ray background 
models predi ct a still-hidd en populatio n of Compton - 
thick AGN ( Treister et all 2004: Wor slev et all [20051 ). 
For obscured AGN (roughly three quarters of the to- 
tal population), the rest frame o ptical light is dom- 
inated by the host galaxy (e.g., iTreister et al.l l2004 
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[Treister fc Urrvl [2OO60 . This obscuring gas and dust 
eventually re-emits the absorbed X-rays at infrared wave- 
lengths. In general the optical/infrared spectral shapes 
of AGN contain two broad components of thermal emis- 
sion: the "Big Blue Bump" seen at optical and ultra- 
violet wavelengths, often attributed to emission from 
a hot accretion disk, visible only in unobscured AGN 
and an infrared excess attributed to re-radiation from 
circumnuclear dust on larger scales visible in all AGN 
in the far-infrared ([Sanders et al.l [19891 : [SandersI [19991 : 
iPier fc Krolig[l99l) . 

Mid-infrared diagnostic studies have found strong con- 
tinuum emission, originating in dust heated by the in- 
tense nuclear ra diation field, a n important indicator of 
AGN emission ([Laurent et al.[ [2000). Indeed, surveys 
with th e Spitzer Space Tel escope Infrared Array Camera 
(IRAC; [Fazio et al. 1 \2004} are finding AGN via color- 
selection techniques by targ eting strong red power-law 
continua dA lonso- Herrero et a l. 20Q6; Hickox et al. 200j; 
[Ponlev et al.[l2007[ ). This emission can be used to find 
even Compton-thick AGN ([Polletta et al.[[2006( ). How- 
ever, many color-selection techniques developed with 
Spitzer data have been based on the infrared colors of 
optically-se lecte d AGN samp les in shallow surveys (e.g., 
Lacv et al.1[2004[ : [Stern et al.|[2005t[Hatziminaoglou et all 



20051) . The ideal AGN selection technique is both highly 



reliable, meaning most objects selected are indeed AGN, 
and highly complete, meaning most AGN are identified. 
To evaluate and improve infrared selection of AGN, it 
is important to investigate how the mid-infrared col- 
ors of X-ray-selected AGN samples compare to those of 
optically-selected samples. 

Finding high redshift AGN in multi-wavelength sur- 
veys requires both a large area and sensitive detec- 
tion limits. The MUltiwavelength Survey by Yale-Chile 
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(MUSYC^; iGawiser et al][200il is a square-degree sur- 
vey of four fields to limiting depths of R ab = 26 and 
Kab = 22 {Kab = 23 in the central areas: [Quadri et aTl 
l2007f ). with extensive follow-up spectroscopy. One 
MUSYC field, the Extended Chandra Deep Field-South 
(ECDF-S), has bee n surve yed by the Chandra X-ray ob- 
servatory ( Giacco ni et al.lT20 02: Al exander et al.l 120031 : 
iLehmer et al.l 120051: IViraniet all 120061) . The sensitiv- 
ity of the Chandra data allows for detection of luminous 
AGN out to large redshifts^ and more modest luminos- 
ity AG N (~ 1042-44 ej-gs s"^) b ack to the quasar era 
[z ~ 2. lBrandt fc Hasingeill2005f l. Recently, the existing 
MUSYC data in the ECDF-S were augmented with very 
deep IRAC Spitzer data in the SIMPLE^ survey (Damen 
et al. 2008, in pre p.) and by the GOODS and SWIRE 
surve ys at 24 ^m (|Dickinson et al.l [20031 iLonsdale et al.l 
120031 ). These two surveys are among the deepest Spitzer 
observations. This results in a large sample of X-ray- 
selected AGN with mid-infrared coverage, unique in the 
combination of area and X-ray and infrared depth. 

In this paper we evaluate mid-infrared selection tech- 
niques using the X-ray selected sample of AGN in the 
ECDF-S. In Section 2, we describe the X-ray and infrared 
data, detail the likelihood ratio-matching technique used 
to find Spitzer counterparts for each X-ray source, and 
present the infrared identifications. We investigate the 
mid- infrared colors of these X-ray sources in Section 3. 
In Section 4 we compare our sources to published IRAC- 
selection criteria and assess the impact of adding the 24 
/im data to the IRAC bands. In Section 5, we discuss 
why the majority of our X-ray-selected AGN would be 
missed by current infrared selection criteria. Conclusions 
are given in Section 6. Throughout this paper we assume 
Ho = 71 km s-i Mpc"\ f2,n = 0.3 and Vlx = 0.7. 

2. OBSERVATIONS & ANALYSIS 

2.1. Chandra Data 

Three X-ray surveys have overlapping coverage of the 
same region of the sky: the very deep central C handra 
Deep Field -South (CDF-S: iGiacconi et all l2002l here- 
after G02; [Alexander et al.l l2003l hereaf ter A03). the 
deep Extended Chan dra Deep Field-Sout h (|Lehmer et al.l 
120051 hereafter L05: IVirani et al.ll2006l hereafter V06), 
and the shallower XMM p ointing that covers this region 
(jStreblvanska et a"ni2004D . For this study, we use only 
the two Chandra surveys because the higher resolution 
of Chandra allows for fainter point source detections. In 
contrast, the XMM data are better suited for X-ray spec- 
troscopy of the bri ghter sources because of the larger 
effect ive area (iStreblvanska et all 120061 : iDwellv fc Pagd 
120061 ). The CDF-S consists of 1 Ms of data reach- 
ing a limiting on-axis flux of 5.2 x 10~^^ ergss"^ cm~^ 
(2.8 X 10-1^ ergss-icm~2) ^j^g 0.5-2.0 keV (2.0-8.0 
keV) band (A03). The 11 individual pointings of the 
CDF-S cover a central area of ~ 300 arcmin^ (G02). The 
four 250-ks Chandra pointings of the ECDF-S overlap 
the deep central region, covering a total area of ^ 1100 
arcmin^. On axis the ECDF-S reaches a limiting flux of 

^ |http : //www . astro ■ yale . edu/MUSYC | 

2 Luminous AGN (Lx ~ 10*" ergs"s-l) can in principle be de- 
tected to very high redshifts in the ECDF-S (z ~ 10), but the 
volume sampled is too small to make such a detection likely, 
http : //www . astro . yale . edu/dokkum/SIMPLE/ 



1.7 X 10"^^ ergscm-^s^i (3.9 x 10"^*^ ergss^^cm^^) 
in the 0.5-2.0 (2.0-8.0) keV band (V06). We combine 
published point source detections from these fields into 
a new unified MUSYC X-ray point source catalog. 

We begin with the most conservative detection cata- 
logs for the CDF-S and ECDF-S, A03 and V06. Since 
the CDF-S goes significantly deeper in flux than the 
ECDF-S and covers a much smaller area between the 
four ECDF-S pointings, we expect that only a fraction 
of sources will be common to both surveys. Both cata- 
logs used similar detection methodologies: first filtering 
the data to eliminate periods of high X-ra y background, 
then using wavdetect ([Freeman et al.l l2002) for source se- 
lection, which is robust in detecting individual sources in 
crowded fields. The MUSYC X-ray catalog first adopts 
the 326 sources presented in A03 Table 3. We next search 
for new sources in V06 Tables 4 and 5 as described in 
what follows. Table 4, the primary catalog in V06 uses 
the same conservative false-posit ive probability th reshold 
for detection (pthrcsh = 10"'') as [Alexander et aP ([2003). 
Given this detection threshold, we expect no more than 
2 or 3 false X-ray sources in either A03 or Table 4 of 
V06. Table 5 of V06, with pthresh = 10-^ allows for a 
larger percentage of false positives but includes a greater 
number of bona-fide X-ray sources in the catalog that 
would otherwise be missed. In order to improve the ab- 
solute astrometric solutions in the C06 catalogs, small 
shifts were applied to the X-ray sources RA and DEC, 
< 1 arcsecond. The shifts were calculated for each of 
the four pointings in the ECDF-S using optical counter- 
parts in the MUSYC BVR image. No shift was necessary 
to ahgn AOS catalog positions. For each source in V06, 
we calculated the distance to the nearest X-ray sources in 
A03 and defined combined X-ray positional uncertainties 

for each pair by cr = \/ '^aos + ^yoe' where (Jxr is the 
1 a error of a Gaussian based on the positional errors 
reported in each catalog. We included in our MUSYC 
X-ray catalog the 503 sources from V06 that were more 
than 4(7 (~ 3 arcseconds) from any source in A03. The 
choice of 4(j is conservative and was determined by de- 
tailed inspections of both the X-ray and MUSYC BVR 
images. The median separation between the reported 
source positions for the 148 X-ray sources found in both 
A03 and V06 was 0.58 arcseconds. A total of 829 unique 
X-ray sources resulted from the union of A03 and V06, 
the first 771 of which were detected with the conservative 
false-positive probability threshold. 

We next investigated the other published X-ra y cata- 
logs in this field for the same two data sets. Lehm er et al.l 
([20051 ). created a larger source catalog for the ECDF- 
S region by retaining data obtained during small X-ray 
background flares and using the same lower significance 
threshold parameter as Table 5 in V06. A very small < 1 
arcsecond positional shift was applied to the L05 cata- 
log to align the astrometry to the MUSYC BVR image. 
Using the X-ray positional uncertainties , we found 141 
additi onal X-ray sources in L05. Finally. [Giacconi et al.l 
([2002D used a set of unique detec tion meth ods including a 
modified version of Sextractor (iB ertin fc Arnouts 199^ 
and running wavdetect with different selection criteria. 
Three catalogs are presented in G02: Table 2 containing 
sources common to both detection methods, Table 3 con- 
taining Sextractor sources only and Table 4 containing 
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wavdetect sources only. A slightly larger fixed positional 
shift (~ 1 arcsecond) was necessary to align the X-ray 
sources presented in G02 with our MUSYC BVR image. 
Using X-ray positional uncertainties, we found a total of 
47 additional X-ray sources not previously presented in 
AOS, V06 or L05: 17 from Table 2, 19 from Table 3 and 
10 from Table 4. Although there is a good possibility 
that some of these additional X-ray sources may be spu- 
rious, those will fail to match to an infrared counterpart 
and thus will not affect the conclusions of this paper (see 
§ 2.3). 

The complete MUSYC X-ray catalog contains 1017 
individual sources whose X-ray properties are given in 
columns 1 through 5 in Table 1. Column 1 presents a 
new unique MUSYC X-ray ID number for each source. 
Column 2 presents the original X-ray catalog and catalog 
ID for each X-ray source (e.g., V06:323 is source number 
323 from Virani et al. 2006). All X-ray source properties, 
such as flux or count rate, are taken from the referenced 
detection catalog. Columns 3 and 4 present the RA and 
DEC, J2000, for each X-ray source, based on the MUSYC 
BVR image. They are on a common system incorporat- 
ing the small ARA and ADEC applied to each catalog 
to shift the coordinates to agree with the MUSYC BVR 
image. Column 5 gives the hardness ratio, defined by 
^^g , where H and S are the X-ray counts in the hard 
(2-8 keV) and soft (0.5-2 keV) bands, respectively. For 
sources detected only in the hard band (soft band), the 
hardness ratio is set equal to 1 (-1). For sources detected 
only in the full band (0.5-8 keV), the hardness ratio is 
set equal to 10. G02 defined the hard band as 2-10 keV. 
Therefore, we adjusted the published hard band counts 
from G02 to equal those that would have been detected 
at 2-8 keV, the hard band defined by the other three 
authors, using PIMMS (version 3.9c; |Mukai|[l993) with 
the average X-ray spectrum measured by G02 (F = 1.4; 
H2-8 = 0.9896 X H2-w). 

2.2. Spitzer Data 

Considerable Spitzer time has been invested in the 
ECDF-S with both IRAC and MIPS (the Muhi-Band 
Imaging Photometer for gpiteer: [Rieke et ai| [2004). The 
GOODS Survey along with additional Guaranteed Time 
Observer (GTO) time gathered deep I RAC an d MIPS 
data on the central region of this field (jDickinson et al.l 
|2003() . Additionally, the Spi tzer Wide-area InfraR ed Ex- 
tragalactic survey (SWIRE; iLonsdale et an[2003l) covers 
a nearly 8 square degree region around the original Chan- 
dra Deep Field-South in a shallow survey mode with both 
the MIPS and IRAC cameras. Finally, a public deep 
IRAC survey was conducted in the region covering the 
ECDF-S (SIMPLE; PI. P. van Dokkum). 

The IRAC data cover 3 — 8 to 5 u limiting depths 
of 24.4 (3.6 pLm), 24.0 (4.5 ^m), 22.2 (5.8 ^im) and 22.2 
(8.0 /im) in AB magnitude. In units of flux density, these 
values are 0.63 /iJy, 0.91 /xJy, 4.7 /iJy, and 4.7 /LtJy, re- 
spectively. A catalog of more than 60,000 sources was 
created with SExtractor using a weighted sum of the 3.6 
and 4.5 //m data (Damen et al. 2008, in prep.). Of these, 
we have reliable photometry in all four IRAC bands for 
38,755 sources (3(7), which will be used in this study. 
The astrometry was calibrated using the MUSYC BVR 
detection image {Rab = 27.1); the resulting positional 
accuracy for individual sources is < 0.3 arcseconds (la; 



Damen et al. 2008, in prep.). 

We also use the SWIRE 24 /xm catalog provided i n 
data release 3 (jLonsdale et al.|[2003HSurace et al.ll2005D . 
which cover the entire area of the X-ray data to AB 
magnitude 18 (0.23 mJy), and the GOODS data in the 
central region, which goes to ~ 5 magnitudes deeper 
(jDickinson et al.l [2003? ). We adopt the aperature fluxes 
from the SWIRE catalog which agree with the published 
GOODS fluxes to better than 5%. The positional ac- 
curacy of the 24 micron data is considerably less than 
that of the IRAC data; however, all 24 micron source 
positions were matched to IRAC sources, decreasing the 
positional uncertainties. 

2.3. Catalog Matching 

Following the work of lBrusa et all (|2005ll2007f l. we use 

the Likelihood Ratio (LR), a statistical matching tech- 
nique, to find Spitzer counterparts to the ECDF-S X- 
ray sources. Although the median positional accuracy of 
Chandra is ^ 0."6, off axis the PSF broadens and be- 
comes circularly asymmetric, leading to larger positional 
uncertainties in faint off-axis sources. The positional ac- 
curacy for each detected Chandra source depends on the 
number of counts detected and the distance of the source 
from the aim point. For several X-ray sources, this po- 
sitional uncertainty encompasses multiple potential in- 
frared counterparts, while for others there may be no 
counterpart to the survey detection limits. Due to the 
high source density of the faintest IRAC sources, there is 
a non-negligible random chance of finding a faint IRAC 
object near an unassociated X-ray source. Simulations 
showed that ^ 60% of random positions in the field fall 
within 4" of an IRAC source, the maximum distance for 
which we find an IRAC counterpart for an X-ray source. 
The Likelihood Ratio takes into account the background 
object density and the positional accuracy of both the 
Spitzer and the Chandra data, and returns the proba- 
bility that a given match is the correct one. A nearest 
neighbor approach, restricted to a typical search radius of 
1.5 arcseconds, would miss 109 IRAC counterparts that 
we recover with our technique. Although larger search 
radii can return more possible counterparts, a nearest 
neighbor approach does not indicate which counterparts 
fall inside the search area by chance and it can only 
distinguish among multiple possible candidates by their 
distance from the reported X-ray source position. Our 
method weights the distance using the positional errors 
and adds a prior on the brightness of the prospective 
counterpart. 

The Likelihood Ratio is the ratio between the probabil- 
ity of finding the true infrared counterpart at a given dis- 
tance from a particular X-ray source divided by the prob- 
abili tv of finding a background object at that distance 
(see [Sutherland fcS aunders 119921 : ICiliegi et all [2003). 
When matching the ith X-ray source, to a possible jth 
infrared counterpart, we calculate: 



where q(m) is our empirically determined probability 
distribution for infrared counterparts as a function of 
magnitude; f(r) is the probability distribution function 
of the positional errors in source i and counterpart j 
(modeled as a two dimensional Gaussian defined with 
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^ gaussian — 

\ ^x,i + ^IR i); and n(m) is the surface den- 
sity of background objects as a function of magnitude. 
When calculating the probability distribution function of 
the positional errors, ux.i is unique to each X-ray source 
depending on its off axis angle and brightness, while no 
such field-dependent distinction is made for the infrared 
counterparts. We adopt ajRAC = 0.3" and cr24^im — 0.5" 
(§ 2.2). Both q(m) and n(m) depend on the depth of the 
infrared data. This technique is inherently Bayesian, and 
the resulting likelihood ratio assumes the priors q(m), 
n(m) and f(r) are known. 

When calculating the background counts of infrared 
sources, nfm), we correct for the fact that there are ac- 
tually fewer faint sources near X-ray objects than near 
random positions in the field (Brusa et al. 2007). This 
follows from the high probability of finding a bright in- 
frared counterpart near any given X-ray position and the 
fact that the fainter sources very close to these bright 
counterparts go undetected. To do this we first estimate 
the background around a sample of IRAC objects with 
a magnitude distribution similar to those IRAC objects 
located within 1" of the X-ray source positions; the po- 
sitional uncertainties make it highly probable that IRAC 
objects within this distance are true counterparts. Sub- 
tracting this estimated background magnitude distribu- 
tion from the magnitude distribution of the actual IRAC 
sources near the X-ray source positions provides a better 
estimate of the magnitude distribution of IRAC sources 
associated with X-ray sources. Experimentation showed 
that using this improved estimate was stable to multiple 
iterations. We then estimate the background around a 
sample of IRAC sources with this final magnitude distri- 
bution for objects associated with X-ray source positions; 
this is the n(m) used in our likelihood ratio calculations. 

Next, we count the total number of infrared objects 
lying within a distance r of all X-ray sources, Ntotai{m) . 
Then we multiply the background distribution by the 
area searched to find these sources; the latter is the num- 
ber of X-ray sources, Nx, times the area of each circle 
searched, irr^ . The difference between N total (m) and the 
expected background, n{'m)Nx'!Tr'^ , is our expected dis- 
tribution of real counterparts, Nreai{m)- It is important 
to understand this is a statistical quantity; rather than 
predicting the expected infrared counterpart for any in- 
dividual X-ray source, it calculates the expected magni- 
tude distribution for all Nx sources. With our expected 
distribution of real counterparts we derive the normaliza- 
tion, Q, as the ratio of the total expected number of real 
counterparts, ^ N^-eai (w) , divided by the total number 
of X-ray sources, Nx- Then we compute our expected 
distribution of true infrared counterparts as a function 
of magnitude: 



q{m) = 



Nreal (™) 



(m) 



(2) 



Finally, we compute the likelihood ratio for each infrared 
object located near an X-ray source using Equation 1. 

When we search for infrared counterparts over a larger 
radius we see more possible counterparts, but the large 
distance between these prospective counterparts and the 
X-ray source drives down the likelihood ratio. We must 
then carefully select the best likelihood ratio which can 



distinguish between true matches and false matches. 
Because the presence of multiple prospective counter- 
parts provides additional knowledge to the likelihood 
ratio, we use the reliab i hty ( R) to select counterparts 
(jSutherland fc Saunderd Il992f ). The rehability is the 
probability that a given infrared counterpart is the cor- 
rect infrared counterpart for each X-ray source. For each 
possible infrared counterpart j, we calculate the reliabil- 
ity by: 

R-^^ns^'^i^. (3) 



Eii?. + (1-Q) 

Here, ^ LRi is the sum over all possible infrared coun- 
terparts for a given X-ray source. The (1 — Q) term 
is the probability of there being no infrared counter- 
part. The calculated reliabilities form a bimodal dis- 
tribution, clustering at values above R=0.9 and below 
R=0.1, and dividing cleanly at R=0.6. Therefore, we ac- 
cept all matches with R > 0.6. The likelihood ratio ap- 
proach assumes exact knowledge of the priors, q(m) and 
n(m). Because we estimate these values from the data, 
the likelihood ratio does not reflect uncertainties in our 
estimation of the background. However, with a threshold 
value of R = 0.6, we find from Monte Carlo simulations 
of the background estimation that fewer than 1% of our 
matches are affected by this uncertainty. 

2.4. Spitzer Counterparts 

We find that more X-ray sources in the ECDF-S are 
detected with Spitzer than at any other wavelength. For 
the full MUSYC X-ray catalog of 1017 sources, we find 
921 Spitzer counterparts in the IRAC data: a 90% re- 
covery rate. Of the 921 sources, 908 have 3a detections 
in all four of the IRAC channels. For many of the X- 
ray sources, blending and confusion of IRAC sources, 
sometimes due to proximity to a very bright source, 
are responsible for the lack of IRAC counterparts. Ex- 
cluding regions where confusion or bright stars interfere, 
we recover infrared counterparts for 98% of the high- 
confidence X-ray source detections, i.e. the subset of the 
X-ray sources detected in the ECDF-S with the conser- 
vative false-probability threshold. For the remaining 2% 
of X-ray sources, either an IRAC source is seen on the 
image near the X-ray position, but this IRAC source is 
not above the source detection limit of the IRAC catalog 
(2 cases), or between 2 and 5 IRAC sources are located 
at larger distances from the X-ray position (2-4 arcsec- 
onds) and no single source has a large enough probability 
of being the correct counterpart (14 cases). We find no 
convincing case of an undetected high significance X-ray 
source in the IRAC image. That is, essentially all X-ray 
sources in our catalog are bright enough in the infrared 
to be detected in the deep IRAC data. 

There are 453 unique MIPS counterparts combing the 
results from the 24 micron surveys in this field. We find 
181 X-ray counterparts in the SWIRE catalog and 318 
counterparts in the GOODS catalog from the central re- 
gion. Combining the counterparts found in the MIPS 
and IRAC data, we have 428 sources with complete SEDs 
from 3.0 - 24 microns. 

All matched counterparts and their mid-infrared fluxes 
are presented in Table 1. Columns 7, 8, 9 & 10 give the 
observed flux in the 3.6, 4.5, 5.8 & 8.0 ^m IRAC bands, 
respectively. The flux units are in /itJy, and are followed 
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by the la errors in parentheses. Column 11 gives the 
MIPS 24 fim flux in fj,3y followed by the la error in 
parentheses. The number is preceded by a G if it comes 
from the deeper GOODS data and an S if it is taken from 
the shallower SWIRE survey. 

2.5. Redshifts 

Although the majority of X-ray sources turn out to be 
AGN, it is also possible for deep X-ray surveys to de- 
tect Galactic stars as well as galaxies that are vigorously 
forming stars. Luminosities above 10^^ ergss""'^ sepa- 
rate out sources that are primarily power ed by black 
hole accretion rather than star formation (Persic et al.l 
120041 : [Lira et al. 2002). In order to study intrinsic X-ray 
luminosities and to apply a luminosity cut to exclude 
Galactic stars and starburst galaxies from our AGN sam- 
ple, we require some knowledge of a source's distance or 
redshift. To this end, we are obtaining optical spectra 
of our X-ray detected sources using Magellan IMACS 
and VLT/VIMOS (Treister et al. 2008, in prep). Thus 
far, we have 281 spectra for our X-ray sample, 110 with 
secure redshifts and 113 with spectral type identifica- 
tions: 42 broad-line sources (BLAGN), 47 non-broad 
line luminous (logL^ > 42 ergs s~^) X-ray sources (non- 
BLAGN), 15 galaxies (logLx < 42 ergss"^) showing no 
AGN feature in their spectra, and nine stars (Treister et 
al. 2008, in prep.). The BLAGN are typically referred to 
as Type 1 AGN, while sources showing narrow-lines are 
referred to as Type 2 AGN. However, these spectral sig- 
natures are often hidden in obscured AGN and moderate 
luminosit y AGN with large host galaxies (|Moran et al.l 
[1)02; bcv ergnmi et al.i r2003). Therefore, we define non- 
BLAGN to be all sources showing narrow emission lines 
suggestive of Type 2 AGN emission plus those sources 
with logLx > 42 ergss"^ that show galaxy-like optical 
spectra. To the 110 redshifts we obtained with IMACS, 
we add an additional 173 sources wi th spectroscop i c red- 
shifts publish ed in the literature fSzokolv et al. '20Q4|; 
Groom et all 120011: iLe Fevre et"a l. 2004; Vanzclla et aL| 
2005( ) . We further include spectroscopic identifications 
from ISzokolv et al.l (|2004( l. which adds 33 BLAGN and 
74 non-BLAGN. Another 274 sources have high qual- 
ity photometric redshifts from the COMBO-17 Survey 
IWolf et al. 2004). We note that our sample of AGN 
spans a range of Lx{0-5 — 8 keV) ~ iq42-45 ergss"^, 
typical for deep X-ray surveys, different from the shal- 
lower survey sources which are on average more lumi- 
nous. This makes a total of 557 X-ray sources (60% of 
the counterpart sample) for which we can calculate the 
X-ray luminosity in the observed 0.5-8 keV band: of these 
416, 76%, have Lx(0.5-8keV) > lO'^^ ergss"! and thus 
are bona- fide AGN. 

In Table 1 , we present the X-ray luminosities and spec- 
tral identifications. Column 12 gives the observed 0.5-8 
keV luminosity calculated from the redshift. For sources 
from G02, no 0.5-8 keV flux was reported. Therefore, 
we adjusted the published 2-10 keV counts from G02 
to equal those that would have been detected at 0.5- 
8 keV, using PIMMS (version 3.9c; Mukai 1993) with 
the average X-ray spectrum measured by G02 (F = 1.4; 
Fo.5_8 = 1.13xF2_io). Column 13 gives the source of the 
redshift estimate used to calculate the X-ray luminosity 
as follows: for VLT/F0RS2 Spectroscopy Version 1.0, 1 
for VIMOS VLT Deep Survey (VVDS) Version 1.0, 2 for 
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Fig. 1. — Color-color plot for the MUSYC X-ray sources. Also 
shown are two SED tracks: (1) the mean type 1 quasar SED from 
Richards et al. (2006) red-shifted from z=0 to z=6 (purple dot- 
dashed line, CCW); (2) pure power-law SEDs, oc {thick pur- 
ple line), with stars marking power-law exponents of a = +1.5 to 
a = —1.5 {bottom left to top right, Aa = 0.5). X-ray-selected AGN 
{light blue crosses) have colors consistent with power-law SEDs 
from the bluest to reddest colors shown here. The black lines indi- 
cate the region suggested by Lacy et al. (2004) for selecting AGN, 
which excludes 42% of X-ray-selected AGN {Lx > W^'^ ergs s^^). 

ISzokolv et all (|200l . 3 for lCroom et"al] ([200lh . 4 for red- 
shifts determined from the MUSYC spectrosc opy (Treis- 
ter et al. 2008, in prep.), 5 for COMBO-17 (iWolf et all 
l2004f ). Those counterparts falling into the categories of 
BLAGN and non-BLAGN are also indicated. 

3. THE MID-INFRARED COLORS OF X-RAY 
SOURCES 

From 3.6 to 8.0 ^m X-ray detected sources show a 
wide range of colors. In Figure 1, a color-color plot of 
the infrared counterparts to the MUSYC X-ray catalog, 
sources with blue colors are in the lower left region and 
those with red colors are located in the upper right. A 
linear relation extends from blue to red along a track oc- 
cupied by power-law SED shapes {fi, ~ z/"). Along this 
track, stars indicate a values of +1.5 (lower left) to -1.5 
(upper right) in steps of 0.5. For each source. X-ray lumi- 
nosities are indicated where known: Lx > 10^^ ergs s^^ 
plotted with light blue x symbols, Lx < 10*^ ergs s~^ 
plotted with light green dashes, and those with unknown 
redshifts plotted with black -I- signs. It is clear that the 
X-ray-selected AGN, Lx > 10*^ ergss""'^, extend along 
the full range of power-law shapes. We also explore the 
infrared colors of the X-ray counterparts identified as 
BLAGN (dark blue boxes) and non-BLAGN (red circles). 
For comparison, we plot the mean Type 1 quasar SED 
compiled by Richards et al. (2006), from redshift to 6 
(dashed purple line, CCW). Unsurprisingly the BLAGN 
lie near the Type 1 quasar track and have the reddest 
infrared colors. Still, many AGN have significantly bluer 
colors than the Type 1 track, falling down and to the left 
on Figure 1. In contrast, the sources classified as non- 
BLAGN are roughly evenly split between red {a < 0) 
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Fig. 2. — Color-color plots for all sources in the MUSYC catalog 
(contours indicate 75%, 50%, 25% and 15% of the IRAC sources). 
The non-X-ray detected galaxies and X-ray-selected AGN are not 
well separa ted in this color space, left The region defined by 
ILacv et "aTl ([200^) includes nearly 20% of all IRAC-detected galax- 
ies. Therefore, this selection region is not effective at identifying a 
sample of AGN. right The selection box of I Stern "etall poos) uses a 
different color-color plot of the IRAC bands. Only ~ 40% of X-ray- 
selected AGN {Lx > 10*^ ergss"^) fall in this region. Because it 
eliminates 90% of the general population of MUSYC sources, this 
selection box is more effective at indicating potential AGN, but 
misses a larger fraction of the X-ray selected AGN. X-ray selected 
AGN occupy a far larger region in mid-infrared color-space than 
optically-selected BLAGN. 



and blue {a > 0) colors. Overall, while the optically 
identified BLAGN tend to have the reddest infrared col- 
ors, the X-ray-selected AGN population is not all red in 
color. 

Next, we compare the IRAC colors of X-ray sources 
(908) to the --40,000 IRAC sources showing no X-ray 
emission to the limits of the published X-ray surveys, 
most of which are inactive galaxies. In Figure 2 we dis- 
play two popular IRAC color-color plots; the boxes in the 
right and left panels correspond to two different proposed 
selection techniques (§ 4) and the contours enclose 75%, 
50%, 25% and 15% of the IRAC sources. The left panel 
shows log (S'5.8)/(S'3.6) versus log (S'8)/(S'4.5), where S is 
the observed flux in /j,Jy, and on the right we plot Vega 
magnitudes: [5.8] — [8.0] vs [3.6] — [4.5]. For clarity we 
plot only Lx > 10^^ ergs s" sources and indicate only 
the BLAGN. On both plots, the power-law SED track 
and Type 1 QSO SED track are repeated from Figure 1. 
Although a larger percentage of the AGN have red colors 
than the galaxies, in neither plot are the X-ray selected 
AGN clearly separated from the rest of the galaxies. 

Because the galaxies and AGN appear to lie along the 
locus of power-laws, we fit the four observed IRAC fluxes, 
3.6-8.0 /ini, to power-law spectral energy distributions. 
Within the short baseline of the IRAC camera, nearly all 
of our sources are consistent with power laws, especially 
given the photometric uncertainties. Each source with 
detections in all 4 IRAC channels was fit by oc j/"; 
94% of all sources and 97 % of the X-ray-selected AGN 
had acceptable fits by power laws (x^, P > 0.05). In 
Figure 3, we compare the histogram of the power-law in- 
dices, a, for the X-ray-selected AGN and other sources. 
Both types of sources show a very broad range in a, but 
the AGN tend to be slightly redder than normal galax- 
ies (more negative a). A KS test rejects the hypothesis 
that the two distributions come from the same popula- 
tion (P=0.0014). While only 17% of the galaxy sam- 
ple have a < 0, 37% of Lx > 10^^ ergss"^ sources fit 
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Fig. 3. — Flistograms of power-law index, a {fi, oc u°'), for 
fits to the Spitzer IRAC photometry. The X-ray sources with 
Lx ^ lO**^ ergs s~^ [solid light blue line) and the subset of X-ray 
sources spectroscopically identified as BLAGN [dark blue cross- 
hatched histogram) have been normalized to the total number of 
X-ray sources with good estimates of redshift and good power-law 
fits (396). Also plotted is the normalized histogram for all MUSYC 
sources well fit by power laws {dotted black line, 38988). In total, 
94% of all galaxies and 97% of the Lx > 10^^ ergs s~^ sources are 
consistent with a power-law SED. A KS-test rejects the hypoth- 
esis that these two distributions come from the same population 
(P=0.0014). Although the AGN are slightly redder than the galax- 
ies on average, both the AGN and the galaxies are fit by a broad 
range of power-law slopes. 



a negative power-law. Optically-selected AGN spectra 
show red continuu m emission at in frared wavelengths, 

with a 1 fe.g- lElvis eraI1ll994[ ). This is consistent 

with our finding that 100% of the BLAGN have accept- 
able power-law fits, 72% with a < (Dark blue cross- 
hatched histogram). As a whole, our power-law fits to the 
BLAGN are consistent with featureless negative power- 
law slopes in the mid-infrared. However, we find a broad 
range of a values even for the BLAGN. Similarly, the 
non- X-ray detected galaxies are bluer as a whole than 
the AGN population, even considering that the sources 
with the reddest a in the X-ray-selected AGN histogram 
may be broad line objects. Indeed, although the X-ray- 
selected AGN have on average slightly more negative a 
values than galaxies, over half (59%) of X-ray-selected 
AGN are fit by power laws with positive exponents simi- 
lar to those of the bulk of the galaxy population. There- 
fore, we find that most X-ray-selected AGN do not show 
red continuum emission at observed Spitzer IRAC wave- 
lengths. 

4. SELECTING AGN BY INFRARED COLORS 

Although selecting AGN through X-ray emission is rel- 
atively unbiased, current X-ray surveys may still be miss- 
ing AGN that are heavily obscured by large columns 
of dust and gas. All AGN, whether obscured or un- 
obscured at optical wavelengths, reprocess and re-emit 
the absorbed radiation in the mid- and fa r-infrared 
(|Granato fc Danesd ll9"9llNenkova et al.l l2002f). Because 
AGN are expected to be luminous at infrared wave- 
lengths, several authors ha ve developed infrared selection 
techn iques for AGN (e.g., iLacv et al.1l2004t iStern et all 
[2005h . Here, we examine what fraction of X-ray-selected 
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A GN would b e selected using these techniques. 

iLacv et all ([2004) found that a sample of Sloan Digital 
Sky Survey quasars occupied a distinct region of mid- 
infrared color-color space (Figures 1, 2). Because the 
selection box was defined using luminous Type 1 quasars 
from the Sloan survey, we look at the subset of the X- 
ray sources that have been identified as broad-line AGN. 
Indeed, 96% of these sources fall inside the selected Lacy 
region. Furthermore, all 61 X-ray sources identified pho- 
tometrically as quasars in the Combo- 17 survey also fall 
inside this box. However, only 58% of all X-ray-selected 
AGN fall inside this region. Thus, the percentage of X- 
ray-selected AGN that lie inside this color selection box 
depends on the number of luminous Type 1 objects in 
the sample. That is, this infrared color selection finds 
the luminous unobscured quasars, but misses a high frac- 
tion of the lower luminosity AGN, many of which look 
like optically normal galaxies. The brightest AGN in 
our sample, Lx > lO^'*'^ ergss"^, are also all located 
inside the Lacy selection box. Additionally, 18% of all 
sources (~ 11, 000 galaxies) lie within the Lacy selection 
box (Figure 2). This is 20 times the total number of 
X-ray sources in this region (529). A reliable AGN selec- 
tion criteria should yield a sample of sources containing 
mostly active galaxies. Because so many IRAC-detected 
sources fall inside this selection box, it is not a reliable 
method for se l ecting AGN in deep surveys. 

IStern et al.l (|2005D define a similar selection region 
using a deeper sample of spectroscopically confirmed 
AGN from the Bootes field. In the right p anel in Fig- 
ure 2 we plot X-ray-selected sources on the IStern et al] 
(|2005f ) color-color plot. Again, we find that the major- 
ity (76%) of our spectroscopically identified broad-line 
AGN fall inside this selection box, but only 40% of X- 
ray-selected AGN do. Over 90% of the brightest AGN, 
Lx > lO^**-^ ergs s~^, are also all located inside the Stern 
selection box. In this case fewer galaxies lie inside this 
box (5163, ~ 8%), but they still out nu mber the 340 X - 
ray sources by a factor of ~ 15. The Ste rn et al.l (j2005f ) 
selection criteria is a slightly more reliable way to select 
for candidate AGN, but at the same time misses a larger 
percentage of the total AGN population. 

Although color selection techniques may work well for 
shallow surveys, where a larger percentage of the AGN 
are bright and show broad lines, they are not as reliable 
in deep surveys. These selection boxes correctly high- 
light the region containing luminous broad-line quasars 
but miss most of the AGN population, in particular non- 
BLAGN and low-luminosity AGN. It is possible that 
some of the galaxies in these selection regions hide lu- 
minous AGN. However, looking at Figure 2, there does 
not appear to be a clear region in color-color space that 
selects a complete sample of X-ray detected AGN and 
avoids returning a large portion of the entire galaxy pop- 
ulation. In an effort to find a more reliable selection tech- 
nique, we now add the 24 /im band data to our analysis 
of infrared colors. 

4.1. 24 fJ.m Data 

Adding the 24 /xm data to a given source's SED more 
than doubles the baseline over which we can determine 
the AGN spectral shape. For ~ 25% of our X-ray sources, 
consistent power-law fits are determined from the 3.6-8 
/xm and 3.6-24 /xm data (examples shown in Figure 4, 
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Fig. 4. — X-ray-selected AGN show a variety of spectral shapes 
below 8 fim, but nearly all 24 /im-detected AGN are very red be- 
tween 3.6 and 24 ^m. In the top plots, a power-law fit to 3.6-8 ^tm 
photometry is significantly bluer than a power-law fit including the 
24 fj,m flux. In the lower plots, the same power-law fits the 3.6-8 
fim and the 3.6-24 fim data. Only ~ 25% of the IVIIPS sample 
falls into the latter category. The plots on the left are for sources 
identified spectroscopically as non-BLAGN while the plots on the 
right are BLAGN. Blue lines show power-law fits to 3.6 — 8fim and 
red lines show fits including the 24 fira data point. 



bottom panels). However, for many sources the best-fit 
power law to the IRAC data alone is much bluer than one 
that includes the 24 /im data (examples shown in Figure 
4, top panels). We found a variety of spectroscopically 
identified BLAGN and non-BLAGN with 24 fim inten- 
sities that confirm the power-law slope from the IRAC 
data and also those whose 24 fim fluxes imply a distinct 
slope. In Figure 4, the sources on the left are identified 
as non-BLAGN from optical spectra while the sources on 
the right are BLAGN. When adding the 24 fim flux to the 
flt, the median non-BLAGN slope changes from a value 
of a3.6-8/im — 0.2 to a3.6-24Aim — —1-0, while the median 
BLAGN slope changes from a value of as.e-S/im = —0.3 
to Q:3.6-24/jm — —1-2. The non-BLAGN appear bluer on 
average in their 3.6-8.0 fim colors, while the BLAGN are 
more likely to show a similar ass-Sfim and as. e - 24pm- 
While nearly all X-ray sources show a negative slope in 
their 3.6 — 24/ir7i spectrum, only for a minority (~ 25%) 
of X-ray sources is it a strictly monotonic rise. 

In Figure 5, we directly compare the fitted index 
cts.e-afim with that of Q!3.6-24/im for all 24 fim detected 
sources. In this figure we include all sources whether 
or not their spectral shape is consistent with a power 
law. The X-ray detected sources span a large range of 
o^3.6-s^ml but in their longer wavelength color, Q;3.6-24/jm, 
they are nearly uniformly red. In contrast, all five galac- 
tic stars detected in X-rays and at 24 fim are in the upper 
left corner (Q!3.6-24pm > 1, dark green points). The his- 
togram of Q3.6-24/im peaks below-1, with '^90% of X-ray- 
selected AGN detected at 24 fim having a3.6-24/im < 
(Figure 5, bottom panel). This histogram is quite differ- 
ent than that for the best-fit as.g-s/im (Figure 5, right 
panel). The vast majority of X-ray sources with counter- 
parts at 24 fim have a^s-2Anm < 0: we question if this is 
because all AGN have negative infrared slopes, or if our 
shallow 24 fim coverage simply misses those AGN that 
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Fig. 5. — We compare the power-law slope (/^ oc v°') from the 
IRAC bands, as.e— Sumi to the power-law slope to all 5 bands be- 
tween 3.6 and 24 /im {main panel). We include all sources whether 
or not fi, oc is a good fit to the data in order to study the ap- 
proximate shape of the SED. For comparison we plot the GOODS 
sample of 24 /xm-detected sources {gray points) and indicate those 
that are consistent with a good power-law fit between 3.6 and 8 fira 
{purple points). The distribution of a3.6-24(/m = peaks below -1.0 
(bottom panel), significantly redder than either of the two peaks 
of the distribution of as.e-g^im for all sources (right panel). The 
cross hatched histogram indicates sources lacking 24 fim detections. 
The deep GOODS data detects the 24 fim fiux for sources with a 
large variety of IRAC spectral shapes. The median uncertainty in 
a for the X-ray sources is indicated in the top left corner. In the 
IRAC bands X-ray selected AGN show a wide range of colors, but 
beyond the observed 8 fim band, nearly all X-ray selected AGN fit 
very red spectral shapes, a < —1.0 However, the 24 /im flux, which 
highlights cold dust, does not help separate AGN from non-active 
galaxies. 

would show blue colors. We indicate sources without 
a 24 fim detection in the cross-hatched portion of the 
a3.6-8/im-liistogram (Figure 5, right panel). The deep 
GOODs data in the central portion returns 24 mag- 
nitudes even for sources with very blue spectral shapes 
in the IRAC bands; we find no significant correlation be- 
tween a3,e-8^m and sources lacking 24 fim detections. 
Compared to 54% of all X-ray counterparts, 48% of the 
X-ray counterparts detected at 24 fim have blue spec- 
tral shapes below 8 fim (as.e-s/im > 0). Only 20% of 
the X-ray selected AGN are consistent with a power-law 
SED from 3.6-24 ^m (x^, p > 0.05), but of those that 
fit a power law, nearly all are negative power laws (34 
out of 38 sources). This contrasts with only 43% of the 
X-ray-selected AGN with IRAC power-law SEDs that 
have as.e-s/jm ^ 0- Although they have a large range 
in color below 8 microns, X-ray-selected AGN detected 
at 24 fim are uniformly red beyond 8 fim. For compari- 
son, we add all 24 fim detected galaxies in the GOODS 
survey to Figure 5 (gray points) and indicate those that 
are consistent with a power-law SED between 3.6 and 
24 fim (purple points). Many of these sources may be 
AGN not detected at X-ray wavelengths, but the rest 
contain significant a mounts of dust wh ich is luminous at 
the observed 24 fim (jEgami et al.ll2004f ). The 24 ^m data 
are not sufficient to differentiate the X-ray selected AGN 
from the non- X-ray detected galaxies. 

5. DISCUSSION 

If the AGN contribution to an X-ray source SED is 
a red power law at infrared wavelengths, the infrared 



emission from most AGN in the MUSYC sample must 
also include a significant contribution from the stars and 
dust in the host galaxy. Stellar spectra are declining in 
the infrared, creating a blue continuum for rest wave- 
lengths A > 1.6 fim. Because the peak of the redshift 
distribution for X-ray-selected AGN in the ECDF-S is 
at z < 1, the AGN host galaxies will typically appear 
blue in the observed 3.6 — S.Ofim band. For the low- 
est redshift objects, contributions from cool dust and 
PAH features can also contribute to the observed mid- 
infrared light (|Laurent et al.|[2"000( ). The 24 fim band is 
sensitive to dust emission in galaxies, making both AGN 
and luminous star forming gala.xies e asily dete ctable 
(iTreister et al.l l2006l : lEgami et all |2004| ; [Le Floc'h et all 
[200ir 

Figure 6 shows X-ray luminosity versus the fitted in- 
frared power-law slope (aa.e-S/im; bottom; as-24fim, top) 
for all MUSYC X-ray sources. Although there is a trend 
for brighter X-ray sources to have redder power laws 
below 8 fim, even some of the most luminous X-ray 
sources in the ECDF-S [Lx = lO*^ - 10^^ ergss"!) 
have a3.6-8/jm > 0. Because sources with similar Lx 
fit such a large range of aa.g-SMirn it is likely that vary- 
ing contributions from the host galaxy (stars and cool 
dust) determine the different IRAC slopes. This is in 
agreement with the large number of non-BLAGN that 
lie outside the infrared color selection box defined by 
optically-selected quasars, and with what is seen in mod- 
estly obscured AGN of Seyfert luminosities, where the 
host galaxy starhght can dominate the optical colors 
(jMoran et al.ll2002l: [Cardamone et al.ll2007[ ). In contrast, 
for infrared power-law selected AGN (as.e-s/im < 0), the 
energetics are dominated by the central AGN and thus 
the signatures are often visible in the optical. Long- 
ward of 8 fim, there is no observed trend of spectral 
shape {as-24fj.m) with X-ray luminosity. Our sample 
agrees with other studies that find that AGN selected 
through IRAC power-law emiss ion are t he most lumi- 
nous of those detected in X-rays (jDonlev e t al. 2007) and 
with studies that look at the most luminous hard X-ray- 
detected AGN and find thern to have red mid-IR colors 
(jCeorgantopoulos et"aIll2007D . 

Supporting the idea that the colors below 8 fim are in- 
fluenced by the host galaxy light is the fact that the X-ray 
spectral shape is entirely uncorrelated with the infrared 
spectral slope but correlated with optical AGN signa- 
tures. We plot a3.6-8A"n and Q;3.6-24/jm versus the X-ray 
hardness ratio, {H — S)/{H -t- 5"), where H is the count 
rate in the hard band (2 — 8 keV) and S is the count rate 
in the soft band (0.5 — 2 keV) in Figure 7. The hardness 
ratio is set to -1 or -|-1 if the source is detected in only one 
of the two X-ray bands, and is not shown here for detec- 
tions only in the full band. The subset of sources with 
optical spectroscopy is indicated by dark blue squares 
(BLAGN) and red circles (non-BLAGN) and sources con- 
sistent with power-law SEDs are indicated (purple stars). 
There are fewer sources consistent with power-law spec- 
tral shapes between 3.6 and 24 fim, but there is no trend 
with X-ray hardness ratio. This result may indicate dif- 
ferent that the X-ray and infrared radiation is produced 
in vastly different environments. The X-rays produced 
near the central black hole are sensitive to obscuration 
near the central region, while the infrared reprocessed 
emission at wavelengths shorter than the observed 24 fim 
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Fig. 6. — Top: At wavelengths longer than the observed 8 ^m, 
there is no observed trend of a8-24(jm with X-ray luminosity. The 
slope of the best fit line to the sources with Lx > 10*'^ ergs s~^ is 
0.04 it 0.06 (dashed line). Bottom: More luminous X-ray sources 
have redder infrared power laws in the IRAC bands, «3.6— 8(j"ii 
although in all cases there is a broad range. The slope of the best-fit 
line to the sources with Lx > lO*'^ ergs s~^ is —0.6 it 0.05 {dashed 
line). It is likely that the galaxy influence, which is proportionally 
larger for less luminous AGN, is the origin of this trend. 

depends on the dust distribution not only around the nu- 
cleus but throughout the galaxy. In contrast, the optical 
identifications based on emission lines are somewhat cor- 
related with hardness ratio: the left side of Figure 7, i.e., 
negative X-ray hardness ratios (relatively soft spectra), 
is preferentially occupied by BLAGN. This reinforces the 
idea that the observed infrared emission is insensitive to 
the geometry, or more likely, includes a substantial con- 
tribution from the host galaxy itself. 

Our results are consistent with the findings for other 
X-ray-selected samples, where a significant fraction of 
the X-ray detected AGN have bluer colors indistinguish- 
able from those of ga la xies (e.g.,[F ranccschini ct al. 2005; 
iBarmbv et all 120061) . iTreister et all f|2006i) concluded 
that a fair fraction of the mid-infrared light in low- 
luminosity AGN arises from the host galaxies, based 
on the large overlap in the infrared luminosity distribu- 
tions of normal galax ies and lower lumin osity AGN in 
the GOODS sample. iHickox et ^ (|2007 ) showed that 
sources were moved into mid-infrared selection boxes by 
combining an AGN template with increasing strength to 
a galaxy template. Mid-infrared selection techniques can 
miss a significant fraction of AGN if the high column 
densities significantly dim the central power-law emis- 
sion and/or star light significantly dilutes the AGN con- 
tinuum. 

Mid-infrared selection techniques can uncover sources 
whose SED is dominated by the central AGN emission. 
AGN undetected in Chandra or XMM surveys because 
of heavy obsc uration are very likel y to be bright in- 
frared sources. lAlonso-Herrero et all ( 2006) used the 24 
/xm-detected sample of sources in GOODS-South to se- 
lect IRAC power-law galaxies. Only half of their sam- 
ple were indi v iduall y detected at X-ray wavelengths. 
iPolletta et al.l ()2006[ ) effectively used infrared-color se- 



FlG. 7. — Power-law indices 0:3.6— 24fim arid aa.g—g^m (top and 
bottom panels, respectively) versus X-ray hardness ratio, {H — 
S)/{H + S), where H represents the counts in the hard Chandra 
band (2-8 keV) and S represents the counts in the soft band (0.5- 
2 keV). Sources that are well fit by a power law over the given 
wavelength range (purple *) are indicated, other sources (black +) 
are plotted with a best-fit value of spectral slopes. The sequence 
of points at -1.0 (-tl.O) are from sources detected only in the soft 
(hard) band. There is no correlation between observed slopes at 
X-ray and infrared wavelengths. The X-ray emission is most likely 
produced close to the black hole, while the mid-infrared emission 
likely comes from dust around the nucleus and in the host galaxy. 

lection techniques to uncover a sample of heavily ob- 
scured quasars. Their selection was based on feature- 
less spectral energy distribution, and they deliberately 
removed any source that could be fit acceptably by a 
galaxy spectral te mplate. A complementary study by 
iDaddi et al.1 ()2007| ) examined non- X-ray-detected galax- 
ies, which were not dominated by power-law continuum 
emission in the IRAC bands, and uncovered a detectable 
contribution from hot dust near an AGN in 20-30% of 
the sample. This emission was discovered in excess mid- 
infrared flux above and beyond the amount predicted 
from star formation indicators at other wavelengths. 
These galaxies lack any detectable signatures of AGN 
emission at optical wavelengths, but stacking results at 
X-ra y wavelengths poin ted towards obscured AGN emis- 
sion (jPaddi et al.ll2007D . 

A well-defined subset of the X-ray counterparts de- 
tected at all infrared wavelengths are well-fit by red (neg- 
ative) power laws from 3.6 to 24 /itm. These sources have 
similar red power laws (a < 0), lie inside the standard 
infrared selection boxes for optically selected AGN, and 
show preferentially broad lines in their optical spectra. 
About 25% of the whole sample fits into this category. 
This is an upper bound since the shallow 24 /im cov- 
erage in the SWIRE survey may bias us towards red- 
der X-ray sources. To check this, we selected a sample 
of sources in the deeper GOODS area that are bright 
enough at 3.6 fim to be detected at 24 /im, even if they 
have a blue slope {a — +1). Of these 26 X-ray sources 
with bright 3.6 /im magnitudes (7713. g < 21.4), 4 have a 
monotomic red slope between 3.6 and 24 /zm. This is a 
slightly lower fraction (15%) but consistent within the 
statistics with the 25% seen in the large sample. These 
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mid-infared power-law sources, tend to be the most lumi- 
nous X-ray sources. Together with their Compton-thick 
counterparts undetected at X-ray wavelengths, these are 
the best candidates for infrared selection techniques. 



6. CONCLUSIONS & SUMMARY 

AGN are luminous infrared sources and deep infrared 
data from Spitzer potentially provide powerful new meth- 
ods for selecting AGN samples. In this paper, we exam- 
ined the infrared properties of a large sample of AGN at 
moderately high redshifts selected through their X-ray 
emission. In the observed 3.6 - 8.0 /im band over half of 
these X-ray-selected AGN show infrared colors consistent 
with those of galaxies. Across the narrow wavelength 
range of the IRAC filters, most galaxies are also consis- 
tent with power-law spectral shapes. As a result, IRAC 
data alone are not sufhcient for a complete or reliable 
selection of (moderate luminosity) AGN. At 24 /zm, dust 
from both star forming galaxies and AGN is detected. 
Unfortunately, observed IRAC to 24 iim colors do not 
effectively separate galaxies powered by AGN from those 
experiencing star formation. Previous studies with IRAS 
have shown the far-infrared to be v aluable for discrimi- 
natin g active from inactive galaxies (jSanders fc Mirabel 
Il996l ): to apply similar techniques to present deep sur- 
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veys requires longer wavelength data, not just at 24 /im 
but at 70 fim and if possible 160 /im, as well. 

To develop complete samples of AGN, we need to un- 
derstand what kinds of AGN each selection criterion re- 
turns. Infrared selection returns a slightly different pop- 
ulation than X-ray selection. Future studies using com- 
plete SEDs from the ultraviolet through infrared will 
help disentangle the contribution from the central AGN 
from that of the host galaxy. We have imaged the entire 
EGDF-S with medium-band filters, which will result in 
high quality photometric redshifts for most of the AGN 
and will allow us to disentangle the nuclear and host 
galaxy light. These data will also reveal the origin of the 
variety of infrared colors shown below 8 /xm. The large 
range in color shown by AGN below 8 potentially 
holds valuable information about the star formation pro- 
cesses going on in the AGN host galaxies and may help 
us understand the AGN-galaxy connection. 
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TABLE 1 

X-RAY AND Infrared Properties of MUSYC ECDF-S X-ray Sources 



ID Cat ID" RA DEC HR,'' S-Qfj-m" S.Ofj.m'^ 24.0/im'' Lx" Notes'" 
(J2000) (A'Jy) (P^y) (orgss-i) 



1 


A03:l 


52.934208 


-27.823917 


-0.09 


1.52 


(0.14) 


1.89 


(0.15) 


1.62 (0.62) 


2.31 


(0.57) 






2 


A03:2 


52.936333 


-27.865639 


-0.43 


6.69 


(0.12) 


7.94 


(0.14) 


5.33 (0.56) 


5.98 


(0.57) 






3 


A03:3 


52.936875 


-27.860889 


-0.24 


11.95 


(0.13) 


14.31 


(0.15) 


17.28 (0.56) 


14.77 


(0.58) 






4 


A03:4 


52.947083 


-27.887028 


-0.34 


15.15 


(0.13) 


22.73 


(0.14) 


32.26 (0.58) 


46.54 


(0.54) 






5 


A03:5 


52.949917 


-27.845972 


-0.54 


5.48 


(0.14) 


4.25 


(0.16) 


3.93 (0.61) 


2.71 


(0.62) 


■ ■ ■ ( ■ ■ ■ ) 43.53 


2 BL 


6 


A03:6 


52.950083 


-27.800528 


-0.47 


23.38 


(0.13) 


28.93 


(0.15) 


30.18 (0.62) 


22.41 


(0.60) 


S: 560.34 (28.32) 41.37 


5 


7 


A03:7 


52.955958 


-27.776222 


-0.37 


13.32 


(0.13) 


16.04 


(0.16) 


20.63 (0.59) 


25.57 


(0.63) 


... ( ... ) 44.17 


5 


8 


A03:8 


52.956208 


-27.842778 


-0.28 


5.33 


(0.14) 


7.03 


(0.16) 


7.80 (0.61) 


7.44 


(0.62) 






9 


A03:9 


52.959958 


-27.845028 


0.10 


15.81 


(0.14) 


19.24 


(0.15) 


15.07 (0.61) 


10.65 


(0.62) 






10 


A03:10 


52.960125 


-27.864389 


-0.42 


12.34 


(0.12) 


11.83 


(0.15) 


8.86 (0.57) 


7.60 


(0.58) 






11 


A03:ll 


52.961000 


-27.883583 


-0.33 


2.42 


(0.13) 


2.85 


(0.15) 


2.54 (0.59) 


6.24 


(0.59) 






12 


A03:12 


52.963167 


-27.847667 


-0.29 


25.68 


(0.14) 


16.94 


(0.16) 


13.30 (0.60) 


8.44 


(0.62) 


■■■ ( ■■■ ) 42.57 


2 iion-BL 



Note. — This tabic is published in its entirety in the electronic edition of the Journal. A portion is shown here for guidance regarding its form and content. 
^ Original detection catalog for X-ray source: and X-ray source ID thcrc-in. Counts and Flux are taken from this catalog.*^ Hardness ratio, defined by , where H and S are the X-ray counts 

in the hard (2-8 keV) and soft (0.5-2 keV) bands, respectively. For sources detected only in the hard band (soft band), the hardness ratio is set equal to 1 (-1). For sources detected only in the 
full band (0.5-8 keV), the hardness ratio is set equal to 10. G02 defined the hard band as 2-10 keV. Theref ore, we adju sted the published hard band counts from G02 to equal those that would 
have been detected at 2-8 keV, the hard band defined by the other three authors, using PIMMS (version 3.9c: 'M ukailll993f) with the average X-ray spectrum measured by 002.^^ Infrared fiux in fijy 
followed by la error in parentheses, from SIMPLE Survey (Damen et al., in prep)*^ 24 micron fiux in fijy followed by la error in parentheses. Source of 24 micron flux is indicated by S for SWIRE 
jLonsdale et al1|2003|) or G for GOODS jDickinson et ani2003l) ^ Observed luminosity [ergs/s] in the 0.5-8 keV band, using references given in fin al column.^ Source of redshift used to calculate the 
lumi nosity: for VLT/F ORS2 Spectroscopy Version 1.0, 1 for VIMOS VLT Deep Survey (VVDS) Version 1.0, 2 for spectroscopic redshifts from fSzokolv et al.| j2004D , 3 f or spectroscopic redshifts 
from[C room et al.| j200lD , 4 for spectroscopic redshifts determined from MUSYC spectroscopy (Treister et al. 2008, in prep.), 5 for photometric redshifts from COMBO-17 jWolf et al.|[200i') . Where 
known, spectra showing broad emission lines (BL) are distinguished from those without broad lines (non-BL). 



